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Abstract 

In this article, we prove the existence and uniqueness of a certain distribution function on 
the unit interval. This distribution appears in Brent's model of the analysis of the binary gcd 
algorithm. The existence and uniqueness of such a function has been conjectured by Richard 
Brent in his original paper pp. Donald Knuth also supposes its existence in |S] where develop- 
ments of its properties lead to very good estimates in relation with the algorithm. We settle here 
the question of existence, giving a basis to these results, and study the relationship between this 
limiting function and the binary Euclidean operator £?2, proving rigorously that its derivative is 
a fixed point of B2. 
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1 Introduction 

If u and v are positive integers, their greatest common divisor (gcd), written gcd(n, v) in the 
sequel, is the largest integer that divides them both. This integer can be computed efficiently 
using a method discovered more than 2200 years ago: Euclid's algorithm. Quoting D. Knuth 
this algorithm is the "grand-daddy" of all algorithms, because it is the oldest nontrivial algorithm 
that has survived to the present day. It is however not the best way to find greatest common 
divisors when dealing with modern computers. In fact, another algorithm, the so-called binary gcd 
algorithm, created by J. Stein (SJ, requires no division but only subtractions, parity testings and 
halving of even numbers (which correspond to shifts in binary notations). These procedures are 
essentially free when compared to the computational cost of divisions. 

The idea of the binary gcd algorithm is basically the following: given two positive integers u 
and v, if halving both numbers is possible at most k time, do it, keeping the values of u and v 
updated. Then repeat the following procedure until both number are equal, say to I: subtract the 
smaller from the greater and when the result is even, divide it by the largest power of 2 possible. 
The gcd of u and v is then I ■ 2 k . This repeated loop will be referred as a "subtract-and-shift cycle" 
in the sequel. 

*The author was supported in part by SNF grant 200021-103683. 



The behavior of the binary gcd algorithm is interesting in several ways. On the one hand, it is 
always important to know the worst case and average case of an algorithm, just from a practical 
point of view, and this is even more important when the considered algorithm has such a wide 
application. On the second hand, the machinery elaborated in order to understand the average 
behavior of the algorithm has led to a deep understanding of it, giving answers as well as rising 
new questions. 

In our case, the worst case the binary gcd algorithm may have to face is a total number of 
subtractions equal to 1 + [log 2 max(u, v)\ ), see, e.g., 

The exact determination of the average behavior of the binary gcd algorithm is however much 
more complex than the analysis of its worst case scenario. Two models have been proposed in 
order to study and analyze the expected behavior of the algorithm. We first describe the model 
created by R. Brent and gives a short description of the model created by B.Vallee at the end of 
this introduction. 

The first model was created by Richard Brent in 1976 pp. In his work, R. Brent exhibits a 
dynamical system describing the binary Euclidean algorithm and provides an heuristic proof of the 
analysis of the algorithm. This dynamical system is described by the binary Euclidean operator B2, 
see (|4.1j) below, that transforms the density associated to the algorithm, step-by-step. However, 
the operator B2 is difficult to analyze, and the question of convergence was left as a conjecture. 
This approach also suffers from the fact that it lies on an unproven connection between a discrete 
and a continuous model, see [Q for more details concerning this last point and 2 for a description 
of the situation 25 years later. 

We now describe this model. Suppose that both u and v, with u > v, are odd, which is the 
case after each subtract-and-shift cycle. Every subtract-and-shift cycle forms u — v and shifts this 
quantity right until obtaining an odd number v! that replace u. Under random conditions, one would 
expect to have u' = (u — v)/2 m with probability 2 _m . This is the heart of Brent's hypothesis. In his 
model, we suppose that u and v are essentially random, except that they are odd and their ratio v/u 
has a certain probability distribution. Let g n be the probability that min(u, v) / max(u, v) is greater 
or equal to x after n subtraction-and-shift cycles have been performed under this assumption. Then 
the sequence of functions {g n }neN satisfies the following recurrence relation 



In the sequel, we will denote by F n {h) the partial sums of the above series. Note that the 
operator F is linear, bounded, since H-FHoo < 2, and that the series converges uniformly for any 
h in C([0, 1]). Computational experiments led Richard Brent to conjecture that the functions g n 
converge uniformly to a limiting distribution g^. Under this conjecture, the function g^ satisfies 
the equality 



g (x) = 1 - x , g n+ i(x) = F(g n ){x) 



where, for all h e C([0,1]) 






and provides the following estimate. If 
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then the expected number of subtract-and-shift cycles in the binary gcd algorithm with starting 
values u and v is \n{uv)/b. 

In this article, we settle this conjecture by proving that the sequence {g n }nm converges uni- 
formly towards a function g^. In order to do so, we first prove that every element of the se- 
quence is convex and differentiable over ]0, 1]. Then we exhibit a compact set of the Banach space 
(C([0, 1]), 1 1.| |oo) which contains g\ and which is left invariant by the operator F defined by (jl.lj) . 
This fact assures the existence of accumulation points of the sequence {<?n}neN, and therefore proves 
the existence of at least one fixed point of the operator F. We prove the uniqueness via an ar- 
gument based on the sequence of derivatives {g' n } n en- On the way, we study the behavior of this 
sequence with respect to the binary Euclidean operator B2 , proving that the sequence converges 
to the unique fixed point of B2 in the L -norm. 

The present work provides a proof that the dynamical system studied by Brent possesses indeed 
a unique limiting distribution. However, it does not shed a new light on the validation of the 
continuous model. In other words, it makes legitimate the work of R. Brent on the analysis of 
the binary gcd algorithm, without validating his model. It also answers a 47-points question of 
D.Knuth [5, p. 355, question 32], who grades the problems of 5 on a "logarithmic" scale from to 



The second work we were referring to is due to Brigitte Vallee who brings another look to the 
situation and leads to a successful analysis using rigorous "dynamical" methods. These methods 
are also the basis for the analysis of several others algorithms [3|S]. In her work, B. Vallee studies 
the operator V2 which describes a slightly different dynamical system. The operator V2 transforms 
the density associated to the algorithm where all the subtract-and-shift cycles are gathered together 
as long as the sign of u — v is constant. As a consequence, the operator V2 is easier to analyze. 
B. Vallee shows that the operator V2 possesses a unique fixed point in some Hardy space and presents 
a spectral gap. She also proves, based on this spectral gap and with the help of a Tauberian theorem, 
the connection between the discrete and the continuous model. Quoting D.Knuth, "her methods are 
sufficiently different that they are not yet known to predict the same behavior as Brent's heuristic 
model. Thus the problem of analyzing the binary gcd algorithm [...] continues to lead to ever more 
tantalizing questions of higher mathematics" . 

Not surprisingly, there is a connection between the two operators B2 and V-z- We refer the 
interested reader to |2j for further details regarding this connection. 

We will use the notation ||.||oo and ||.||i for the supremum norm and the L 1 -norm of functions 
defined over [0,1], and log 2 for the logarithm in base 2. Let us recall that a series of function 
^ n>0 /i n (x) verifies the so-called Weierstrass criterion (see |3J III. 4]) over a subset A of R if we 
have 



Let us mention that the notation for our g n and g^ are different in both and 0- Brent 
uses F n and F^ and Knuth uses G n and G. 

2 Convexity and Regularity 

We prove in this section that the elements of the sequence {g n }n€N are convex and decreasing 
functions over [0, 1] and differentiable over ]0, 1]. 



50. 




sup |/i n (x)| < oo- 

x<=A 



(1.2) 



n>0 
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Let m be a 2 x 2 matrix with real coefficients a, b, c and d . Such a matrix acts naturally on 

via 

a b 
c d 



mix 



. . ax + b 
(x) 



cx + d 



and this action satisfies (mi • m,2)(x) = m\{rri2{x)) for all pairs m\,m2 of 2 x 2 matrices, where 
• is the usual matrix product. From now on, we will identify a matrix m with the real function 
associated to it. Let us define the set M as follows: 

M = < m : [0, 1] -> R | m(x) = with a, 6 > 0, c, d > , ad - be / 

Note that for any element m of M, sgn (m) := signum(ad — 6c) is well-defined, since a common 
factor at the denominator and the numerator of m does not affect the sign of ad — be. This function 
satisfies the equality sgn (mi o mz) = sgn (mi) • sgn(m2) for all pairs mi,m,2 in M. Let us also 
define the set S as the set of all series 2~}ieN e * m *> where m; G M, satisfying the following three 
points: 

1. £j = ±1 and Si ■ sgn (m«) < , Vi G N. 

2. The series verifies the Weierstrass criterion over [0, 1], i.e., 

'^2e i m i £S =^> ||mj||oo < +oo. (2.1) 

3. The following series converges: 

k_ 

Cidi e%x di 



E\didi — biCi\ , . OjX + 6j 
: < +oo , where mAx) = -. (2.2) 



Note that the series ()2.2j) is well-defined since a common factor at the denominator and the 
numerator of rrn does not affect the terms of the series. A typical element g of S can be written as 

^ = a^rh (23) 

Cj 3^ "T" Ct? 

For sake of clarity, let us recall two facts about series of functions: first, if a series of functions 
satisfies the Weierstrass criterion (|1.2[) on a set, then it converges uniformly and absolutely on it, 
and the limit does not depend on any permutation of the sum. This result applies also for double 
sums. Second, if the derivatives of the partial sums of a convergent series of function converges 
uniformly then the series is differentiable and its derivative is the limit of the derivatives of the 
partial sums, see, e.g., Thm. 2.13, Thm. 2.9, Thm. 4.3 and Thm. 6.18 of [I]. 

The definition of the set S takes its roots in the following two lemmas, which are the keystones 
of the article. 

Lemma 2.1 Every function g in S is a convex, decreasing, continuous function over [0, 1] and 
continuously differentiable over any compact set o/]0, 1]. 

Proof: A function m in M is convex and decreasing if and only if sgn (m) < 0. Indeed, we have 

'ax + b\' ad -be n , , 

< <S=^ ad-bc< 0, 



cx + d J (cx + d) 2 
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and 



ax + b 
cx + d 



-2c- ad ~ b " >0 



ad — be < 0. 



(cx + d) 3 

Using the first two points of the above definition, any element g in S is a uniform limit of convex, 
decreasing and continuous functions, and is therefore convex, decreasing and continuous. Let us 
prove now that any element of S is continuously differentiable over any compact interval of ]0, 1]. 
Let < e < 1. For g as in (|2.H|) . the partial sums of g' satisfy 



' N 

E 



aiX + bi 



N 



\ a%di b{C{ 
^ {ax + di) 2 



\i=0 V i=0 

and the definition of S shows that this series satisfies the Weierstrass criterion over [e, 1] since for 
all x G [e, 1], 

-2 I O'idi b{Ci | 



< — — - < e 



(cix + di 



(cj£ + di) 2 



(Ci + di) 2 



and 



yields 



1 

4cd 



i , ,x 2 < — Vc, d > 0, 
(c + d)^ 



biCi\ ^ e 



i>0 



(cj + dif 



'z ^ v — *\ \(Xjdj — bjCj\ 

— -E - 

i>0 



Thus, the partial sums of derivative converge uniformly over [e, 1] to a limiting function which is 
the derivative of g. This finishes the proof. 



Lemma 2.2 Let F n be the partial sums of the series hl.l)) . If g : [0, 1] 

then F n {g) G S for allnGN and F{g) G S. 



Proof: Let us define the following particular elements of M: 



Hk{x) 





2 k 



(x) 



and Vk(x) 





2 k 



(x) 



□ 



is a function in S, 



1 



2 k x + 1' 



x + 2 fc 

Note that these functions map the interval [0, 1] in itself and satisfy 

sgn(/Xfc) > and sgn(z^) < 0. 

Let us prove that if <7(x) = J2ien £iTni{x) is in S 1 , then F(g) lies inside S. The proof for the partial 
sums F n (g) is similar, although infinite sums might become finite. We have 



F(g)(x) 



Y,t~ k (aMx)) - g{v k (x))) 

k>l 

E 2 ^ ^timiiVkix)) - y^gjmt(i/ fc (a;)) J 

k>l VieN iGN / 



(2.4) 



Based on the Weierstrass criterion (|2.1j) . we have 

iin^iioo < E 2 ~MEn 

° Mfcl loo 



fc>l \«GN 

s eWe 

fe>l VieN 



+ ^ ||m, o VkWoo 



m, ; 



+E 



rn, 



(2.5) 
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and therefore the double sums in Q2.4|) can be rearranged in any simple sum 

F(g){x)=J2 £ i M i( x ) 

where EiMi(x) is either of the type Ej • 2~~ k mj(iJik(x)) of the type — Ej ■ 2~ k m,j(v k (x)). Clearly, 
F(g) has the correct structure to be an element of S. We must now prove that this function fulfills 
the three points of the definition of the set S. Inequality (|2.5f) shows that the latter series fulfills 
the Weierstrass criterion, directly proving the second point. Since 

Ej • sgn (2~ fe • (rrij o /j, k )) = Ej • sgn (rrij) • sgn (/i fc ) = Ej • sgn (rrij) < 0, 

and 

Ej • sgn (-2~ k ■ (rrij o v k )) = -Ej ■ sgn (rrij) ■ sgn (v k ) = Ej • sgn (rrij) < 0, 

the first point is verified. Let us check the validity of the third point. A straightforward computation 
shows that if g is as in (j2.3j) then the analogue series as ()2.2j) for F(g) is the following double series: 

v - ^ 2\<Xidi — biCi\ 
2^2-^ (d + dMi ■ 2 k ' 

k>l i£N v 1 1,1 

This series is convergent since the inequality (cj + d()di > C{d{ yields the following estimate: 

\ ^ \ ^ 2\aidi — biCi\ 1 \&idi — biCi\ 

(c i + d i )d i -2 k < 2^^k'2^ ^ 

k>l ieN v 1 k>l i€N 

\ -< \CLidi — 6jCj 

= 2 z. — ^ — <+0 °- 

This shows that the third point is fulfilled and the lemma is then proven. □ 

Proposition 2.3 Every element of the sequence {Sn}neN\{o} * s i n Thus every element of the 
sequence {g n }neN is a convex, continuous and decreasing function over [0,1], continuously differ- 
entiate over ]0, 1]. 

Proof: Since go(x) = 1 — x, go fulfills the conditions of the claim. The function g\ is as follows 
Si 0*0 = Y, 2 k \i + 2 k x ~ x~+¥j = ^ \2 k + 2 2k x ~ 2 k x + 2 2fc 

k>l v 7 k>\ v 

and a straightforward computation shows that the function g\ above is an element of S. Lemma 
12.21 shows by induction that g n is an element of S and is therefore convex, decreasing, continuous 
over [0, 1] and continuously differentiable over any compact subset of ]0, 1] by Lemma |2~T1 □ 



3 Existence of an accumulation point 

In the current section, we prove that the sequence {g n }n^N possesses at least one accumulation 
point in the Banach space of continuous function defined over [0, 1], with the supremum norm. Let 
us define the following two subsets of this Banach space: 

K x = SH {g\ g(0) = = 0}, (3.1) 

K 2 = {g G C([0, 1]) 1 1 + 3/2 • x log 2 x - 5x < g(x) <l-x}, (3.2) 
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where S is the closure of S in the supremum norm and log 2 is the logarithm in base 2. Note 
that any element of S is a decreasing, convex and continuous function, being a uniform limit of 
such functions. The definition of K 2 seems odd at first sight. The key point is that a function in 
K 2 cannot come close to 1 with a too steep slope when x goes to 0. We start with the following 
proposition: 

Proposition 3.1 The operator F verifies the following properties: 

1. F(Ki) C K u 

2. F(K X n K 2 ) C K 2 , 

and therefore F(K\ Pi K 2 ) C Kif]K 2 . 



Proof: The map F being continuous, we have F(S) C F(S). Lemma \'2.2\ implies that F(S) C S, 
therefore F(S) C S. The fact that F (g)(0) = 1 and F (g)(1) = when g(0) = 1 and g(l) = is 
straightforward. This proves the first point. 

Suppose g is a function in K\C\K 2 . The inequality F(g)(x) < 1 — x is obvious since, F(g) being 
an element of S, is convex and lies below the secant joining (0, 1) to (1, 0). It remains to show that 

F(d)(%) > 1 + 3/2 • xlog 2 x — 5x. 
Based on the definitions of F and K 2 , we have 

F(9)(x) = ^22- k (g(x/(x + 2 k ))-g(l/(l + 2 k x)) 

k>l 

> E 2 ~ fc (i + r • K * - io g2 (x + 2*)) ) 
fc>i 



X , 1 

-5-— T-1 + 



x + 2 k l + 2 k 



x 



^•xW x V 1 1 3 x V 1 lQ g^ + 2 ^ 

2 62 \Z^2 k x + 2 k 2 \^2 k x + 2 k 



'1^2*' x + 2 fc ) + I S 2^ ' l + 2 fc x ) 
\fc>l / \fc>i / 



Note that, for x G [0, 1], we have 

9fc _|_ Ofc — Ak Q 



2 fc x + 2^ - 4^ 
fc>i fe>i 



and 

1 lQg 2 (x + 2 fc ) < ^ 1 = 
Ofc ™ 1 Ok — 2-^i Ok 



2 k x + 2 k ~ ^ 2 k 

k>l k>l 

Based on Mellin's transform, the equality 

E Y k • TTW-x = 1 + x log2 x + x - P(log2 x) + \ " B" 1 )* aJ^ry** ( 3 - 3 ) 

fc>l k>2 
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where 

2-7T ^— ^ sin27rA;y 
^ V > = m~2 ' ^ sinh(2/br 2 /ln2) 

can be proven. A proof can also be found in p. 644], where it appears as the main step in the 
computation of the function g\. As a matter of fact, the function P{y) is small, and can be bounded 
in absolute value by 8 • 10~ 12 , c.f. 5 . We will however only need a far less accurate bound. Since 
sinh(i) > e*/4 for t > In 2/2, we have 

| P(y) | < *1 . y U e -^/^) k = . = 1.5549.. • HT 11 < 1/4. 

1 yyM In 2 ^\ ) In 2 e 2 ^ 2 / lri2 - 4 ' 

k>l 

For x £ [0, 1] , the terms of the alternating sums on the right-hand-side of ()3.3j) decrease in absolute 
value. This sum can therefore be bounded above by its first term, 2x 2 . Since xlog 2 x < over 
[0, 1], the previous estimation of F(g) becomes 

3 13 
F{g)(x) > - • xlog 2 x • - - - • x • 1 

1 1 X 2 

—5 • x ■ — + 1 + x log 2 x — x-- + — — 2x 

3 , 35 2 
= 1 H — ■ x log 9 x x — 2x 

2 52 12 

3 /25 
= 1 + — • x log 2 x — 5x + I — x — 2x' 



3 

> 1 + — • x log 2 x — 5x. 

This last estimate finishes the proof of the proposition. □ 



Note that the proof of the previous proposition also shows that if a function g is convex and 
in K2, then F(g) is in K2 as well. Indeed, the only property needed from K\ in the proof that 
F{K\ n K2) C K2 is the convexity of elements in K%. Let us state this result as a corollary: 

Corollary 3.2 If a function g : [0, 1] — * M is convex and in K2, then F(g) is in K%- 

We turn now to a result of compactness. 
Proposition 3.3 The set K\ n K2 is compact in the Banach space (C([0, 1]), 1 1 ■ 1 1 00 ) - 

Proof: The set K\ is closed in (C( [0, 1] ) , 1 1 - 1 1 00 ) being the intersection of two closed sets. The set 
if 2 is clearly closed as well. Consider the set of Holder functions over [0, 1] (with parameter 1/2). 
These are the functions / for which 

N m \f(x)-f(y)\ , 

N1/2U) = SU P — j [172— < °°- 

x ^y \x y\ 

Based on an argument of equicontinuity, it can be verified that the set 

K A ,B = {/ € C([0,1]) I H/IU < A, N 1/2 (f) < B} 
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is compact in C([0, 1]) for any A, B > 0, see e.g., Chap. 4, Sect. 6 of Let us show that K\ (IK2 C 
if 1,5. The only non-trivial point to be checked is the fact that if g £ K\ n K2, then A r 1/ / 2 (5 f ) < 5. 
The function g being decreasing and convex, we have, for < y < x < 1, 

\g{x) -g(y)\ g{0)-g(x-y) 1 — g(h) 

< = -j= — , with h = x — y > 0. 



[x-y] 1 / 2 ^Jx-y y/h 

Using a property of the elements of K2, we also have 



1-^) < -3/2/, kg, ft + 5* = ^ ^ 



The maximum value of the latter function, defined over [0, 1] , is reached for h = 1 and therefore 

5Vh- 3/2V^log 2 /t 



\g{x) -g{y)\ < 



h=l 



Taking the supremum, we obtain the expected result. The set K\ n K% being a closed subset of a 
compact metric space, it is itself compact. This proves the proposition. □ 

The previous two propositions give directly the next corollary, since any compact set in a metric 
space satisfies the Bolzano- Weierstrass condition: 

Corollary 3.4 The sequence {flVijneN possesses at least one accumulation point in K\ n Ki. 

Proof: The function go is convex and in K2. By Corollary 13 .2\ F(go) = gi is therefore in K2. 
This function is also in K\ (see the proof of Proposition 12. 3(1 and thus any element of the sequence 
{ffnjneN but go is in the compact K\ n K2- The conclusion follows by the Bolzano- Weierstrass 
property. □ 



4 Behavior of the derivatives and uniqueness of the accumulation 
point 

In the current section, we prove that the sequence {g n }nen hi fact possesses only one accumulation 
point goo. This proves that the sequence converges to this well-defined function in K\ n K2 since 
a sequence in a compact metric space with only one accumulation point converges to this point. 
In order to achieve this goal, we study the sequence of derivatives {g' n } n eN in the topology of the 
L 1 -norm over ]0, 1]. Then, based on a property of the binary Euclidean operator B2, defined below 
by (|4.1|l . we show the uniqueness of the accumulation points of both the sequences {g n }neN and 

The sequence of derivatives {<7^}neN does not converges uniformly, or at least we do not know 
it, and therefore nothing tells us that an accumulation point g^ possesses a derivative which is the 
uniform limit of the derivatives of the subsequence. However, g n being convex for all n, we will see 
in the sequel that the sequence of derivatives does converge but in a weaker topology, the topology 
of ^([0,1]). Recall that any convex function h over [0, 1] is differentiable almost everywhere in 
[0, 1] and that it is absolutely continuous pj, i.e., 

h(x) = h(0) + ( h'(t)dt, Vx G [0, 1]. 
Jo 

The following lemma sheds light on the convergence of the derivatives of convex functions: 
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Lemma 4.1 Let {/ n }neN be a sequence of convex functions defined over [a,b], differentiable over 
]a,b[, and converging uniformly to a function f. Then: 

1. If E is the subset of point of]a,b[ where f is differentiable, then for all xq in E, the sequence 
{fn(xo)}neN converges to f'(x ). 

2. If the functions f n are monotone, then the sequence {fn} n eN converges to f in the L l -norm. 

Proof: Let xq E E. There exists ho > such that [xq — ho, xq + ho] C]a, b[. The functions f n being 
convex, we have VneN, VO</i</io, 

f n (x - h) - fn{x ) .,. . fnjxp + h) - fnjxp) 
h ^/nWS h 

When n goes to infinity, this leads to the following inequalities V < h < ho , 

f(x -h)- f(x ) f(x + h)- f(x ) 

< liminf f n {xo) < limsup/ n (xo) < r . 

n n n n 

Taking the limit when h goes to 0, we finally have 

f'(x ) < liminf /^(xo) < limsup/^(x ) < f'{x ), i.e., lim f n {x Q ) = f'(x ). 

n n n—*oo 

In order to prove the second point, note that since / is convex, the set [a, b]\E has measure 0, and 
thus f' n converges to /' almost everywhere. Without loss of generality, suppose the functions f n 
are increasing, i.e., f' n > almost everywhere. The functions f n and / are absolutely continuous, 
and thus 

lim f \f n \= lim f f n = lim (f n (a) - f n (b)) = f(a) - f(b) = f f'= [ \f'\. 

A direct application of the dominated convergence theorem shows that f n converges to /' in L , 
see also Ex.21, p. 57 of 0. □ 



Consider the following linear operator, obtained by taking the formal derivative of the series 



(dU): 



B 2 (h)(x) = J2 



k>l 



+ 2 k 



h 



+ 2 k 



+ 



1 



l + 2 fc ; 



h 



1 



1 + 2 fc ; 



(4.1) 



This operator is referred as the "binary Euclidean operator" in the literature. It was first 
studied by Richard Brent 1 . It is not clear at first sight for what class of function the operator 
Bi should be defined. If we consider its action on L 1 ([0, 1]), then the operator B 2 is a contraction 
with respect to the Li-norm: 

2 



r \B 2 {h)(t)\dt < e tijT 

Jo k>1 Jo V + 



2 k 



h 



t 



+ 



1 



E 2 



— k 



l + 2 k t 

l/(l+2 fc ) 



t + 2 k 
1 



k>l 

r-1 



1 + 2 k t 
\h(y)\dy + 



dt 

dt 
i 

l/(l+2 fc ) 



\Ky)\dy 



\h(y)\dy. 



(4.2) 
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The first equality comes from the changes of variables y = t/(t + 2 k ) in the first integral and 
y = 1/(1 + 2 k t) in the second integral. As a consequence, the operator B 2 can be defined over 
the entire Banach space L 1 ([0, 1]), and this operator is continuous with respect to the topology 
generated by its norm: 

B 2 : L\[0, 1]) — » L\[0, 1]) and B 2 G £(^([0, 1]), L J ([0, 1])). 

This property was already noticed by Brent in his original article PQ. Here is a first application of 
Lemma 14. II 

Proposition 4.2 If h is a function in K\, c.f. iTO) . then F(h)' = B 2 {h') in L 1 ([0, 1]). 
Proof: Let h £ S D {g\ 5(0) = 1,5(1) = 0}. Let us define 




By Lemma 12.21 and 12 . 1 1 these functions are convex, decreasing and continuously differentiable over 
]0, 1]. Therefore, over ]0, 1], we have 




The sequence {/^IneN, being the partials sum of the series that defines B 2 , converges to B 2 (h') 
in L 1 ([0, 1]). The condition of Lemma f4. II are fulfilled and since the sequence {/ n }neN converges 
uniformly towards F(h), we have F(h)' = B 2 {h') in -L 1 ([0, 1]). 

In general, if h £ S n {g\ g(0) = 1,5(1) = 0}, there exists a sequence {/i n }neN of function in 
S H {g\ 5(0) = 1,5(1) = 0} that converges uniformly to h. We use again Lemma 14.11 to have that 
h' n — ► h' in L . Since every h n belongs to S Pi {5)5(0) = 1,5(1) = 0} the partial result above 
applies and thus F(h n )' = B 2 (h' n ) for all n G N. Since {F(h n )}n£N is a sequence of decreasing 
convex functions that converges uniformly to F(h), thanks to the continuity of F with respect to 
1 1 . 1 1 oo , we can once again apply Lemma 14.11 to this sequence. Taking the limit leads to the result 
since F(h n )' — ► F(h)' in L l and B(h' n ) — ► B(h') in L 1 as well, because of the continuity of B 2 
with respect to ||.||i. □ 



Based on these properties, we can prove the following expected theorem, using another time 
Lemma 14. II 

Theorem 4.3 The sequence {5 n } n eN possesses a unique accumulation point, and therefore con- 
verges uniformly to a limiting function g^ which is a fixed point of the linear operator F. The 
sequence of derivatives {g' n } n £N converges almost everywhere and in the L^-norm to g'^, which is 
a fixed point of the linear operator B 2 . 

Proof: Let us consider 500, an accumulation point of the sequence {5 n }ngN- Based on Proposition 
14.21 we have 

g' 00 = F(g OD y = B 2 (g' 00 ) in ^([0, 1]). 
Consider the following sequence of non-negative real numbers: 

Un= [ Woo ~ 9n\dt = \\g'oo ~ g' n \\l , 71 G N. 

Jo 
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Then, using the previous equality, Proposition 14.21 and the fact that B 2 is a contraction, we see 
that the sequence {it n }neN is decreasing: 



Un+l 



< 



a'oc 



g'oo 



9n+l\\l 

F(g n y\\i 

£ 2 (<4) - B 2 {g' n ) 
B 2 {g' 00 -g l n )\\ 1 

/ _ /|| 

9 00 9n 1 1 1 



Let {g nk }k&N be a subsequence of {g n }neN that converges to goo. Note that the conditions of Lemma 
14. II are fulfilled and therefore the sequence {g' nk }k&N converges in L ([0, 1]) and almost everywhere 
to g'oo. This implies that the decreasing sequence {ii n }neN possesses a subsequence that converges 
to and therefore 

lim u n = 0. 

In other words, the sequence {g' n } n eN converges to g'^ in L 1 ([0, 1]). Thus, since 
lim g n (x) = lim I g n {0) + / g' n (t)dt 



1 + / 9'oo(t)dt = goo(x), 




we see that the sequence {g n }n£N converges point-wise to goo. This makes impossible the existence 
of another accumulation point. As explained at the beginning of the section, this shows that the 
sequence {g n }n€N converges to g^- D 



5 Conclusion 

Theorem 14.31 shows that the operator B 2 has a unique eigenfunction with eigenvalue 1. Compu- 
tational experiments show that the next eigenvalues seem to be conjugate complex numbers A 
and A close to 0.1735 ± 0.00884vi, with |Ai| = |A 2 | = 0.1948. Therefore, B 2 seems to present the 
spectral gap Vallee's operator V 2 possesses. The method described in this article does not seem to 
extend in such a way that this spectral gap can be proved. This would imply the exponential speed 
of convergence already suspected in . 

We did not prove that the function g'^ is continuous over ]0, 1], which is strongly suspected. If 
proven, this continuity would directly imply the uniform convergence of the sequence of continuous 
and increasing functions {g' n }n€N over any compact set of ]0, 1] because of a theorem of Dini. 
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